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ABSTRACT 

The dark matter content of the Universe is likely to be a mixture of matter and antimatter, perhaps 
comparable to the measured asymmetric mixture of baryons and antibaryons. During the early stages 
of the Universe, the dark matter particles are produced in a process similar to baryogenesis, and dark 
matter freeze-out depends on the dark matter asymmetry and the annihilation cross section (s-wave 
and p-wave annihilation channels) of particles and ant ip articles. In these rj— parametrised asymmetric 
dark matter models (77ADM), the dark matter particles have an annihilation cross section close to the 
weak interaction cross section, and a value of dark matter asymmetry rj close to the baryon asymmetry 
t]b- Furthermore, we assume that dark matter scattering of baryons, namely, the spin-independent 
scattering cross section, is of the same order as the range of values suggested by several theoretical 
particle physics models used to explain the current unexplained events reported in the DAMA/LIBRA, 
CoGeNT and CRESST experiments. Here, we constrain rj— parametrised asymmetric dark matter by 
investigating the impact of such a type of dark matter on the evolution of the Sun, namely, the flux 
of solar neutrinos and helioseismology. We find that dark matter particles with a mass smaller than 
15 GeV, a spin-independent scattering cross section on baryons of the order of a picobarn, and an 
rj— asymmetry with a value in the interval 10~ 12 — 10~ 10 , would induce a change in solar neutrino 
fluxes in disagreement with current neutrino flux measurements. This result is also confirmed by 
helioseismology data. A natural consequence of this model is suppressed annihilation, thereby reducing 
the tension between indirect and direct dark matter detection experiments, but the model also allows a 
greatly enhanced annihilation cross section. All the cosmological rj— asymmetric dark matter scenarios 
that we discuss have a relic dark matter density flh 2 and baryon asymmetry r/B in agreement with 
the current WMAP measured values, Q. DM h 2 = 0.1109 ± 0.0056 and r\ B = 0.88 x lO" 10 . 
Subject headings: dark matter-elementary particles-stars:evolution-stars:interiors-Sun:interior 
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1. INTRODUCTION 

During the last three decades, overwhelming evidence 
has been found for the existence of dark matter in the 
Universe. This achievement is the result of a careful 
analysis of the available cosmological observational data, 
combined with a variety of well structured theoretical 
physical models coming from quite different and com- 
plementary research fields in particle physics, cosmology 
and astrophysics. Such studies have led the way to identi- 
fication of the basic gravitational effects of dark matter, 
and their contribution to the formation of structure in 
the universe. In spite of this, the fundamental nature of 
the dark matter particles remains a mystery. 

Current observational studies suggest that the matter 
present in the Universe is compos ed predominantly o f 
dark matter particles and barvons (|Munshi et al.ll201ll) . 
Recent measurements of the cosmic microwave back- 
ground by the WMAP team ([Larson et al.l 1201 lh give 
precise measurements of the dark matter density and 
the baryonic density: ilnAih 2 = 0.1109 ± 0.0056 and 
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fl B h 2 = 0.02258± u ; uu gf£. These observational studies 
also measure the imbalance between baryons and an- 
tibaryons, i.e., the baryon asymmetry, which is found 
to be equal to (0.88 ± 0.021) x 10~ 10 . 

Among the most popular candidates for dark 
matter are a group of particles that occur natu- 
rally in supersymmetric extensions of the standard 
model of particle physics, usually called Weakly 
Interacting Massive Particles (WIMPs), the neu- 
tralino, the lightest supersymmetric stable particle 
being the typical example. Significant constraints 
in the properties of WIMP candidate s and similar 
parti cles have been made usi ng starsdBertone et all 
20051). such as th e Sun (e.g.. lLopes fc Silkl l2010allbl; 



Taoso et all l2010t ICuinberbatch et alj I2010D, sun 



stars (Casancllas & Lodcs 201 I tH) and neutron stars 
IKou 



(e.g.. [Kouvarisl l27bTl) iKouvaris fc Tinvakovl l2?)llL [2010; 
iBertone fc Fairbairnl [2008D . Although, WIMP candi- 
dates can be either dirac or majorana particles, usu- 
ally, WIMPs are considered to be of the later type, as 
such these particles do not have an asymmetry such as 
the baryon asymmetr y. Recently, severa l authors (e.g 



Gudnason et~afl l2006t iFoadi et al.1 120091: iHooper et al l 
| 2005h followin g previous work (e.g. iNussinov"! 119851 : 
iKaplan I Il992f) have proposed a new type of matter, 
known as asymmetric dark matter, which has an asym- 
metry identical to baryons ([Kaplan et al.1 120091 ) . Like 
WIMPs, such particles are non-relativistic massive parti- 
cles that interact with baryons on the weak scale, thereby 
having a sizeable sc attering cross-section with baryons 
([Kaplan et al.l |2009T >. but unlike WIMPs, such a type 
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of dark matter is produced in the primordial Universe 
by a mechanism similar to baryogenesis. This type of 
dark matter, much the same as baryons, is considered 
to have an asymmetry that we choose to represent by 
the parameter rj. Hence, this dark matter is composed 
of an ungual amount o f matter and antimatter (e.g. 
iFarrar fc~Z ahariias 2006). Furthermore, th ese particles 
have a mass of the order of a few GeV (e.g. iKang et all 
1201 It ICohen et all 12010( 1 . We will refer to this type of 
dark matter as 77— parametrised asymmetric dark matter 
(7/ADM) or simply as rj— asymmetric dark matter. 

Several experiments committed to direct dark matter 
searches show evidence of positive particle detection, al- 
though these results are still very controve rsial and not 
universally accepted: the DAMA/LIBRA (iBelhni et al. 
2012tlBernabei et al.ll2008h and CoGeNT (|Aalseth et al. 
2011(1 experiments find evidence of a particle candidate 
with a mass of the order of a few GeV (likely between 5 
and 12 GeV), and a spin-independent scattering cross 
section of baryons of the order of 10~ 40 cm 2 . The 
CRESS (H rown et ail l2012f) experiment also points to 
unexplained events consistent with direct detection of a 
light mass particle. Unfortunat ely, other direct dete ction 
experiments, suc h as CDMS (lAhmed et al.l 1201 lh and 
XENON10/100 (lAprile et al.l 120 lit I Angle et al.l boilh 
find no indication of the existence of such particles. Nev- 
ertheless, several theoretical explanations based on the 
existence of asymmetric dark matter have been suggested 
to explain and acco mmodate all these positive and nega- 
tive detections (e.g.lChang et al.ll201CtlFarina et al1l2011l : 
iHooper fe Kelsoll2011t iDel Nobile et al.ll2011ll . ~ 

In this paper, we investigate the origin of such light 
77— asymmetric dark matter particles, and discuss how 
this type of matter influences the evolution of stars. 
We look for their impact on the present structure of 
the Sun. The Sun, by means of two groups of observ- 
ables, helioseismology and solar neutrino fluxes, pro- 
vides one of the most powerful t ests of stellar evolution 
(jTurck-Chieze fc Couvid at 2011) and o f alternative the- 
ories of modern physics and cosmology ([C asanella s et all 
120121 iTurck-Chieze et al.ll20Ll Hopes et al l l2002t l . Al- 
though both groups of observables are equally relevant 
for such types of studies, we have chosen preferentially to 
study the impact on solar neutrino fluxes, as neutrinos 
are the most sensitive probe of changes in the structure 
of the solar core. In particular, we will compare the neu- 
trino fluxes of a "Sun" evolving in a halo of asymmet- 
ric dark matter with the current measurements of solar 
neutrino fluxes, namely, the a B and 7 Be neutrino fluxes. 
Moreover, we will complete the study with a succinct 
helioseismology diagnostic. 

In the reminder of this paper, we consider that the 
dark matter asymmetry is identical to the baryon asym- 
metry, and likewise leads to an unbalanced amount of 
particles and antiparticles. This rj— asymmetry occurs 
well before the epoch of thermal decoupling of the dark 
matter. We do not discuss here any of the possible mech- 
anisms for the generation of such asymmetry, but rather 
treat this quantity as a free parameter. This asymmetric 
dark matter framework is used to study the impact of 
such matter in the evolution of the Sun for a wide range 
of particle masses, annihilation cross sections and dark 
matter asymmetries. 

In Sec. 2 we present the basic properties of 



rj— parametrised asymmetric dark matter and explain 
how the relic dark matter density of the present day 
Universe is computed. In Sec.(s) 3, 4 and 5 we dis- 
cuss the changes caused in the Sun by the presence of 
77— parametrised asymmetric dark matter in its core, as 
well as the impact on the flux of solar neutrinos and 
helioseismology. In the last section, we summarize our 
results and draw some conclusions. 

2. ASYMMETRIC DARK MATTER IN THE EARLY 
UNIVERSE 

In the early stages of the Universe, the cosmic fluid 
rapidly reached thermal equilibrium, as a consequence of 
standard and dark matter particles having a very high 
rate of interactions. In this work, we allow dark matter 
to be constituted by a mixture of particles and antipar- 
ticles, x an( i X with a mass m x , and q v being the num - 
ber of internal degrees of freedom (jDrees et al.l 120061) . 
The Universe is considered to have a temperature T and 
an effective number of relativistic degrees of freedom g+ 
(|Iminnivaz et al.l 1201 lh . Unlike symmetric dark matter, 
X and x particles are not necessarily identical. Following 
the usual convention, the asymmetric nature of the dark 
matter is defined by the parameter 7/ which is equal to the 
difference of particle and antiparticle populations. With- 
out loss of generality, we conventionally define rj > 0, i.e., 
particles are more abundant than anti-particles. Here 
rj stays constant throughout the evolution of the Uni- 
verse. The asymmetry rj is similar to the asymmetry for 
baryons rjs, both of which originate in the early phases 
of the Universe. If not stated otherwise, in the numerical 
examples, the population of x an d X particles is consti- 
tuted by elementary particles with mass m x ~ 10 GeV, 
g+ = 90 and g x — 2. This is equiva lent to a particle 
population forme d by Dirac fermions ((Pent et al.l 120101: 
iDrees et al.ll2006l) . 

As the Universe expands, the interactions between all 
particles become sparse, the temperature of the plasma 
drops, the Universe gets cooler, up to the moment that 
the XX annihilation rate drops below the Hubble expan- 
sion rate, leading to dark matter becoming decoupled 
from the rest of the cosmic fluid. The dark matter den- 
sity, FIdm > has frozen out and has been constant ever 
since. The relic densities of x an d X particles are de- 
termined by solving the coupled system of Boltzmann 
equations that follows the time evolution of the number 
density o f particles and antipa rticles in the expanding 
universe (jlminnivaz et al.ll201lT ). 

The present-day dark matter density, FIdm, depends 
on the relic densities of particles and antiparticles, Q x 
and Q x , which in turn depend on the mass of the particle, 
m x , the XX annihilation cross section, a and the dark 
matter asymmetry rj. 

Primarily, the value of Udm relies on the properties 
of the XX annihilation cross section. In most cases, the 
thermal averaging of the XX annihilation cross section 
times the relative velocity of colliding particles v can be 
expanded as (crv) = a + bv 2 + C(w 4 ). If the annihilation 
initial state is unsuppressed, the first term of the previ- 
ous equation dominates (a ^ 0). This process is known 
as the s-wave annihilation channel. Alternatively, if the 
initial state is suppressed, a = and b ^ 0, this is the 
p-wave annihilation channel. This expansion is valid for 
all known examples of s-wave channel, p-wave channel or 
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Fig. 1. — Asymmetric dark matter particle models: The figure shows iso-curves of the relic dark matter density Q^Af^ 2 as a function 
of the asymmetry parameter r) and the annihilation cross section (av): (a) s-wave annihilation channel (a ^ 0,6 = 0) and (b) p-wave 
annihilation channel (a = 0,6 7^ 0). The results are for dark matter particles with m x = 10 GeV, g x = 2 and g t = 90 (see text). In the 
figure, indicated with blue lines are the r\— asymmetric dark matter particle cosmological models compatible with the present day dark 
matter density. In the figure, indicated with blue lines are the asymmetric dark matter particle cosmological models compatible with the 
present day dark matter density. This cosmological model is computed for two types of annihilation channels, (a,rj) or (b,rj) (see text). 
We consider the observational window to be such that 0.10 < ^Idm^ 2 < 0.12 (or —1.0 < log^oM^ 2 ) < —0.92). This corresponds to the 
region between the blue lines. The observed value of Q.nMh 2 is shown as a green line. The current measurements of the cosmological density 
and baryonic asymmetry are il DM h 2 = 0.1109±0.0056 and r) B = (0.88±0.021) X 10" 10 , or log(tt DM h 2 ) ~ -0.9551 and log(?y s ) ~ -10.06 
ULarson et al.ll20Tll; IHutsi et afll201lT) . 



both channels, up to an accuracy of a few percent. The 
expression (av) can be simplified further, if the reacting 
particles are non-relativistic, (av) = a + Qbx^ 1 + 0{x ) 
where x defines the ratio of ro Y over t he temperature of 
the Universe T (|Iminnivaz et al.ll2011f) . 

In the case of symmetric dark matter particles (x = x)j 
the relic abundance il x is mainly determined by the an- 
nihilation rate (crv). Furthermore, fi x is equal to the 
dark matter value Qdm- Alternatively, in the case of 
asymmetric dark matter (x 7^ x): X and X particle con- 
tributions have to be added such that Q dm — &x + j 
as there are more particles than antiparticles (or the re- 
verse). The antiparticles are annihilated away more effi- 
ciently, with large numbers of particles left behind with- 
out a partner to annihilate them. Consequently, the relic 
dark matter abundance is determined not only by the 
(av) as in the symmetric dark matter case, but also by 
the asymmetry parameter 77. 

With the objective of computing the evolution of the 
Sun in different cosmological scenarios, we start by deter- 
mining the relic densities of particles and antiparticles, 
f2 x and f2 x for specific 77— asymmetric dark matter mod- 
els. The computation is done by following c losely the 
numerical procedure of llminnivaz et all (|2011[ ) . Figure 1 
shows the dark matter density VLdm for light asymmet- 
ric particles with m x ~ 10 GeV computed for several 
values of 77 and (av). Two types of XX annihilation rates 
are considered, pure s-wave and a pure p-wave annihila- 
tion channels. Annihilation channels with non- vanishing 
a and b terms are qualitatively similar to the previous 
ones. 

In general, such results show that there are two limit- 
ing asymmetric dark matter scenarios, one for very high 
and other for very low values of rj: (i) for the high values 
of 7/, the dark matter is "strongly" asymmetric as the 
relic density is dominated by particles over antiparticles, 
il x il x w £Idm', this scenario is identical to the baryo- 
genesis process related to ordinary baryons. (ii) for the 



low values of 77, the dark matter is "weakly" asymmetric 



or symmetric (77 « 0), Q x 



Oojvf; the final relic 



abundances of particles and antiparticles are compara- 
ble. This generally corresponds to the standard thermal 
WIMP scenario. 

If we consider the present measurements of dm 
(jHiitsi et al.l[20Tll) . we find that the number of asym- 
metric dark matter scenarios that are compatible with 
observations is relatively reduced, i.e., only 77— (av) dark 
matter models with a relic density 0.1053 < Udm^ 2 < 
0.1165. Nevertheless, for reasons of convenience 8 , in this 
work we choose the observational interval to be such that 
0.10 < SlDMh 2 < 0.12 (region between blue lines in Fig- 
ure 1). 

If we restrict our analysis to the asymmetric mod- 
els of particles with a mass of 10 GeV compatible with 
dark matter density observations, it is possible to de- 
termine a critical value of (av) c for which there is a 
change of dark matter regime (see figure 1). The crit- 
ical value (av) c is approximately 1.6 x 10~ 24 cm 3 s _1 and 
1.2 x 10 -23 cto 3 s _1 for s-wave and p-wave annihilation 
channels. There is also a maximum value of asymme- 
try ?7 C , above which the asymmetric dark matter mod- 
els have a £Idm larger than the current observational 
value. ?7 C is equal to 3.6 x 10 -11 for the s-wave annihi- 
lation channel and 2 x 10~ n for the p-wave annihilation 
channel. If a model has (av) < (av) c , Q.dm becomes 
uniquely dependent on 77 (independent of (av) ). This 
corresponds to cosmological models for which the XX an_ 
nihilation rate is so efficient that the relic abundance 
depends only on the initial 77 asymmetry. Conversely, 
if a model has (av) > (av) c , the dark matter density 
becomes purely determined by the value of (av) (inde- 
pendent of 77). It follows that the dark matter density 
content is determined by the low value of the annihila- 

8 The recent measurements of QrjM^ 2 are such that 0.1053 < 
Cl DM h 2 < 0.1165 or -0.9776 < log(n DM h 2 ) < -0.9337, never- 
theless for clarity of argument we choose a slightly larger interval 
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tion rate. 

It is worth noticing that the dependence of SI dm with 
the particle mass is quite small for asymmetric dark mat- 
ter models of light particles (m x < 20 GeV). In cosmo- 
logical dark matter models for particles with a mass be- 
tween 5 GeV and 20 GeV, (av) c for both channels is iden- 
tical to the case of a particle with a mass of 10 GeV. For 
the same range of masses, rj c varies between 8 — 2 x 10~ n 
for a s-wave (a ^ 0) channel and 7 — 2 x 10~ n for p-wave 
channel (b 0). 

3. CAPTURE OF 77— PARAMETRISED ASYMMETRIC 
DARK MATTER BY THE SUN 

The accumulation of particles (x) and antiparti- 
cles (x) inside the star depends on the star's grav- 
itational field and the interaction of the dark mat- 
ter particles with baryons. In our description of 
the impact of the rj— asymmetric dark matter on 
stellar evolution, we follow closely earlier discussions 
of the impact of "classical" asymmetric dark mat- 
ter 9 with the evolution of the Sun and other stars 
(IScott et alJl2009HCumberbatch et"aIll2010HTaoso et al.1 
120101 ICasanellas fc Lope s 2009; lLopes et al. 1 120111) . 

In the following we discuss the capture and accumula- 
tion of rj— parametrised asymmetric dark matter by the 
Sun. Although, such particles interact with baryons in 
a similar way as WIMPs, the fact that we have two dis- 
tinct populations of particles, leads to a quite different 
impact on the evolution of the star. Like WIMPs, the 
amount of rj— asymmetric dark matter captured by the 
star depends explicitly on the mass of the dark mat- 
ter particle m x , the cross-section for scattering with 
baryons, namely, the spin-dependent scattering cross- 
section asp , and the spin-ind ependent scattering cross- 
section as 1 (jLopes et al.ll2011l ). However, unlike WIMPs, 
rj— asymmetric dark matter, depends also on the two 
channels (p-wave and s-wave) for annihilation (av) and 
the 77— asymmetry parameter. The total number of par- 
ticles N x and antiparticles N x that accumulates inside 
the Sun at a certain epoch is computed by solving the 
system of coupled equations: 



dNj 
dt 



a - C a N x N x 



(1) 



with i being \ or X- The constant Cj gives the rate of 
capture of particles (antiparticles) from the dark matter 
halo, and C a gives the annihilation rate of particles and 
antiparticles in the Sun. This approach is entirely dif- 
ferent from other studies of accumulation of dark matter 
inside stars, in part because previous work only studied 
the capture of the population N x of dark matter parti- 
cles. In such cases, the system of coupled equations (TTJ 
is substituted by one single equation. In the following, 
we will consider dark matter particles with a mass larger 
than 5 GeV, for wh ich the evaporation of particles is neg- 
ligible (|Gouldl lT990') . Furthermore, we neglect the rate of 
capture of dark matter particles by scattering off other 

9 In the literature, we find another description of asymmetric 
dark matter, that is distinct of the type of dark matter discussed in 
this work. Although both types of asymmetric dark matter accu- 
mulate by identical processes inside the star, their effects in the evo- 
lution of stars is rather different. We decided to call r\— asymmetric 
dark matter to highlight the difference 



dark matter pa rticles that ha ve already been captured 
within the Sun (|Zentnerll2009| ). 

The capture rate for particles and antiparticles by 
the Sun's gravitational field and their scatterings off 
baryons is computed num erically from the integral ex- 
pression of iGouldl (l!987f) . implemented as indicated in 
iGondolo et al.l ( 2004 ). The description of how this cap- 
ture process is imp lemented in our code is discussed in 
lLopes et al.1 (|2011l ). The scattering of particles and an- 
tiparticles on baryon nuclei are identical. It follows that, 
similarly to the symmetric dark matter case, as d is only 
relevant for hydrogen nuclei; and asi is important for 
the scattering of dark matter particles on heavy nuclei. 
If the value of a si is larger or equal to as d , the capture 
of dark matter particles is dominated by the collisions 
with heavy nuclei, rath er than by collisions with hydro- 
gen (jLopes et al.ll201lD . 

If not stated otherwise, we assume th at the density of 
dark matter in the halo is 0.38 GeVcm~ 3 (|Catena fc Ulliol 
l2010f K the amount of particles and antiparticles in the 
halo is in the same proportion as in the local Universe, 
the stellar velocity of the Sun is 220 kms -1 and the 
Maxwellian veloc ity dispersion of dark matter particles is 
270 kms" 1 (e.g.. iBertone et al.|[2005h . Recent measure- 
ments of local dark matter density yet to be published 
propo se a value of 0.3 GeVcm" 3 (jBoyv fc Tremaine I 
mm and 0.85 GeVcm" 3 (jGarbari et al 1120121 ). The an- 
nihilation rate C a is computed from 



C a = N~ N x J (av){r) n x (r)n x (r) 4ttV dr, (2) 

where (av)(r) is the \X annihilation rate as defined in 
the previous section, and n x (r) and n x (r) are the number 
density of particles and antiparticles. 

In this paper, in contrast to previous work, we assume 
that the dark matter particles and antiparticles present 
inside the Sun annihilate by the same physical process as 
the one that occurs in the early Universe, therefore, i.e., 
(av) — a + Qbx^ 1 + (D(x~ 2 ) where x is now the ratio of 
mass of the dark matter particle over the local tempera- 
ture of the Sun's plasma. It follows that inside the Sun 
the s-wave (a ^ 0) annihilation channel is identical to the 
one present in the early Universe, but the p-wave (b ^ 0) 
annihilation channel is quite different, as the tempera- 
ture inside the Sun is larger than the temperature in the 
early Universe. 

As the star evolves, after some time, t is larger than 
t~a = y C a C x , the rate of accumulation of particles in- 
side the Sun N x is proportional to the difference between 
the capture rate of the particles and antiparticles, i.e., 



N x = D x with D x = C x 



C x . It follows that the total 



number of particles and antiparticles is given by N v ^ 
iVv + D x t and N x » C x /{C a D x t) dGriest fc Seckell 
1987). The value of ta depends of the values of dark 
matter parameters. For a typical dark matter halo of 
density 0.38 GeVcm -3 constituted by particles with a 
mass of lOGeV and annihilation rate 10 _24 cm 3 s~ 1 , ta of 
the order of 1 million year. Particles and antiparticles, 
once captured by the Sun during its evolution, end up 
in thermal equilibrium with baryons. The particle (and 
antiparticle) population follows a Maxwellian velocity 
distribution. The number density of particles (antipar- 
ticles), rij(r) N ie - m x<P(r)/T a wnere y c in the central 
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Fig. 2. — Asymmetric dark matter model impact on the solar neutrino fluxes: The asymmetric dark matter particles have a mass 
of 10 GeV and a scattering cross section of baryons such that ergo = 10 -40 cm 2 and osi = 10 — 36 cm 2 . The asymmetric dark matter 
cosmological models have a value of (av) and r) for which the cosmological models have Qjjj\,[h 2 value consistent with the observational 
window, i.e., 0.10 < CIdm^ 2 < 0.12 or —1.0 < log (Hjjm'i 2 ) < —0.92(08). This corresponds to asymmetric dark matter models for which 
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r/i 2 is in between the blue lines Figure 1 (see text). 



temper ature of the star and <p(r) i s the gravitational po- 
tentia l (fCasancllas & Lopes 200^; iGiraud-Heraud et al.l 
H990h . 

4. IMPACT OF ^-PARAMETRISED ASYMMETRIC DARK 
MATTER ON THE SUN'S EVOLUTION 

The code used to compute the evolution of 
asymmetric dark matter is a modifi ed ver s ion o f 
the stellar evolution code CESAM (jMorell I1997D . 
The basic reference model without dark matter 
is calibrated to produce a solar standard model 
(|Turck-Chieze fe Lopesl 11993 ) identical to others in the 
literature (iGuzik fc MussacMl201Cl: iSerenelli et al.l 120091: 
iBahcall et al.ll2005bHTurck-Chieze et al.ll2010D . The mi- 
croscopic physics (updated equation of state, opacities, 
nuclear reactions rates, and an accurate treatment of mi- 
croscopic diffusion of heavy elements) are in full agree- 
m ent with the s t andar d picture. The solar mixture 
of lAsplund et al.l (|2005| ) is used in the computation of 
models. The modified evolution models are calibrated 
to t he present solar radius, luminos ity, mass, and age 
(e.g. iTurck-Chieze fe Couvidatl 1201 ll ). The models are 
required to have a fixed value of the photospheric ratio 
(Z/X) & , where X and Z are the mass fraction of hydro- 
gen and the mass fraction of elements is heavier than 
helium. The value of (Z/X)^ is determ ined according 
to the solar mixture lAsplund et al.l (|2005l ) . 

The star evolves from the beginning of the pre-main 
sequence until its present age. Each solar model has 



more than 2000 layers, and it takes more than 80 time 
steps to arrive to the present age. For each set of dark 
matter parameters, a solar-like calibrated model is ob- 
tained by automatically adjusting the helium abundance 
and the convection mixing length parameter until the 
total luminosity and the solar radius are within 10~ 5 
of the present solar values. Typically a calibrated solar 
model is obtained after a sequence of 10 intermediate 
solar models, although the models with a large concen- 
tration of dark matter need more than 20 intermediate 
models. The increase of the number of iterations is re- 
lated with the rapid variation of the structure in the 
center of the Sun. The values of the calibrating parame- 
ters, mixing-length parameter and the initial content of 
Helium changes slightly relatively to the standard solar 
model, depending upon the accumulation of dark mat- 
ter in the core of the star. The maximum variation ob- 
tained for these parameters due to the accumulation of 
dark matter in the centre of the Sun is an increase of the 
mixing-length parameter of 5% and the initial content of 
Helium is decreased by 0.5%. 

The r\— asymmetric dark matter particles impact the 
evolution of the Sun, by a physical process identical 
to WIMPs. Likewise they provide an effective mecha- 
nism for the transport of energy, for which the efficiency 
depends locally on the average mean free path of the 
dark matter particles between successive collisions with 
baryons. Two distinct regimes of the transport of en- 
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Fig. 3. — Comparison of the sound speed radial profile between the solar standard model and different solar models evolved within 
r\— asy mmetric dark matter halos: The red dotted curve corresponds to the difference between inverted profiles and our solar standard 
model (Basu ct al. 2009; Turck-Chiczc ct al. 1997). The dark matter particles have a mass of 10 GeV and a scattering cross section o 
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ergy are usually considered: local and non-local trans- 
port corresponding to small and large average mean free 
paths. Both energy transport mechan isms are imple- 
mented in our code (jLopes et al.ll20lil ). This physical 
pro cess leads to the reduction of the temp erature gradi- 
ent (jLopes et al.ll2002tlLopes fc Silkll2010bD . In the more 
extreme case, the high frequency of collisions between 
bary ons and dark matter particles forms an isothermal 
core (|Lopes fc Silkll2n0l . 

This new type of asymmetric dark matter can have a 
much larger effect on the evol ution of the Sun tha n the 
usual symmetric dark matter (jLopes fc Silldl2010bl ). be- 
cause the star can accumulate a much larger amount of 
particles (and antiparticles) in its core than in the latter 
case. This is due to the fact that asymmetric dark matter 
depends, among other parameters, on two major param- 
eters that determine the concentration of dark matter 
inside the star: dark matter asymmetry, which deter- 
mines the unbalanced amount of particles and antiparti- 
cles; and the annihilation cross-section, which establishes 
the annihilation efficiency of particles and antiparticles. 
This is consistent with the results shown in Figure 1. 

Figure 2 shows the variation for the 8 B and 7 Be neu- 
trino fluxes relative to the solar standard models. The 
large reduction in S B and 7 Be neutrino fluxes is a direct 
consequence of the reduction of the temperature inside 
the solar core. 

This diminution of temperature is also visible in the 
sound speed profile near the centre of the Sun. Fig- 
ure 3 shows the profile of sound speed computed for a 
few rj— asymmetric dark matter scenarios. It also shows 
the observed sound speed profile obtain from the helio- 
seismolog y data of the BISO N and GONG observational 
networks (Basu et al. 2009). This sound speed is con- 
sistent with a previous sound speed inversion computed 
from high accuracy data obtained by the GOLF and MDI 
instru ments of the SOHO mission (|Turck-Chieze et"aTI 
1997). In the solar core, the impact of rj— parametrised 



dark matter on the sound speed is much larger than can 
be accommodated by the curren t solar standard model 
(|Turck-Chieze fc CouvidatJ 120111 ). Moreover, although 
the sound speed is not inverted i n the very central re- 
gion of the Sun (jBasu et al.ll2009| ). the reduction in the 
sun's speed profile caused by the presence of dark mat- 
ter is much larger than previo usly tentative sound speed 
inver sions of the Sun's core (jTurck-Chieze fc Couvidati 
1201 lh . As shown in figure 3, these dark matter mod- 
els produce a decrease in the square of the sound speed 
of the order of 3%-5%, well above 1% of the uncer- 
tainty of the current solar model. This reduction in 
temperature is accompanied by significant changes in 
the Sun's core structure, which leads to visible effects 
on other type of seismic d iagnostics, such as the small 
acoustic mode separationdLqpes fc Turck-Ghie^ H99l 
iTurck-Chieze fc Lopeslll993D . ICumberbatch etaIT(|2010[ ) 
have found that dark matter particles that accumulate 
inside the Sun's core and produces a temperature varia- 
tion of this order of magnitude, also have a small acous- 
tic mode separation quite different from helioseismology 
data. 

The relationship between the solar neutrino fluxes and 
the sound speed with the temperature is easily estab- 
lished, if we notice that the 8 B and 7 Be neutrino fluxes 
depend on the central t emperature T c as T c 24 and T c 10 
(Ba hcall fc Uhnen[l99l iBahcalJ I2002T) . and the square 
of sound speed can be expressed as Cl ~ T c / n c where 
ji c m olecular weight in the centre (jTurck-Chieze fc Loped 
1993). This different sensitivity to temperature explains 
why the diagnostics provided by &( 8 B), $( 7 Be) and 
are quite distinct. We have also computed the solar neu- 
trino fluxes of other pp-nuclear reactions and the results 
follow a similar behaviour. Other seismic diagnostic are 
possible, such as comparing the observed acoustic oscil- 
lations to the theoretical acoustic oscillations, or their re- 
spective small separations. Nevertheless, the sensitivity 
of seismic quantities to changes in the Sun's core struc- 
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ture is always less visible than in solar neutrino fluxes. 
5. DISCUSSION 

The current st andar d solar model 

(ITurck-Chieze & Loped 119931 : iTurck-Chieze et al.l 

l2010t ISerenelli et al.l 1 2 llh is in agreement with the 
$( 8 _B) and $( 7 _Be) neutrino fluxes measur e d by 
SNO and Borexi n o de t ector (lAharmim et al.l 120101 : 
iBellini et al.1 [2fTTfl l20lil 120121 : lArpesella et alT f2C)08h . 
The current S B neutrino flux observational deter- 
mination in the case of no-neutrino oscillations (or 
electron neutrinos; including the theoretical uncertainty 
of neutrino flux solar neutrinos in the solar stan- 
dard model) is $( 8 B) = 5.0 5^;^ x 10 6 c m ~ 2 s~ 1 
for the SNO experiment (Aharmim etaf] [2010) 
and <f>( s B) = 5.88 ± 0.65 x 10 6 cm - 2 s~ 1 for the 
Borexino experiment ((Bellini et al.l 12010). The Borex- 
ino experiment measures 7 Be solar neutrino flux 
$( 7 Be) = 4.87 ± 0.24 x 10 9 cm~ 2 s~ 1 , under the as- 
sumption of the MSW-LMA scenario of solar neutr ino 
oscillations (|Bellini et alj|2011t lArpesella et al.ll2008D . 

The $( 8 _B) neutrino measurements made by Borex- 
ino and SNO experiments suggest that the Sun's core is 
slightly hotter than expected, since the $( 8 _B) measured 
value is higher than the value predicted by the current 
solar standard model. This discrepancy is validated by 
the Borexino measurement of $( 7 Be) which is sensitive 
to a region slightly off the Sun's centre. The first mea- 
surement of the pep neutrinos was done recently by the 
Borexino experimen t, $(pep) = 1.6 ± 0.3 x 10 8 cm _2 s _1 
(|Bellini et al.l 120111 ). The experimental value also sug- 
gests that the Sun is hotter than the solar standard model 
predicts. This diagnostic is quite reliable once we realise 
that §(pep) is strongly dependent on the luminosity of 
the star. All these independent solar neutrino experi- 
ment results suggest that the solar standard model is 
cooler than the actual Sun. 

The internal structure of the Sun is well known 
by means of solar neutrinos and helioseismology data. 
Therefore, the theoretical uncertainty of the standard 
solar model is consistently taken into account by com- 
paring the predictions of these two probes with obser- 
vations. Several authors have shown that an impor- 
tant source of uncertainty on the calculation of the so- 
lar neutrino fluxes (electronic flavour) comes from the 
unclear measurem ents of heavy elemen t abundances in 
the S un's surface (|Bahcall et al.l 2005b. a: lAsplund et al.l 
2009). One other source of uncertainty is related with 
the determinat ion of several pp-reaction rates and elec- 
tron s creening (Mussack & Dappen 20111 : ISerenelli et al.l 
1201 It IWeiss et al.ll2001l ). Similarly, other standard and 
no-standard physical processes that contribute for the 
evolution of the star can also be an important source of 
uncertainty. Among others we can refer to the follow- 
ing indirect processes: convective overshoot, low-z ac- 
creti on, mass loss, solar rotation and meridional circula- 
tion dGuzik fc Mussack 120101: ITurck-Chieze et "all 120101: 
ITurck-Chieze fc Couvidadl2011l ). These suggestions have 
been made mostly to resolve the discrepancy between the 
theoretical s ound speed as obtained wit h the new solar 
abundances (jAsplund et al.1 12009L l2005| ) and the sound 
speed obtained from the acoustic oscillations. As pre- 
viously mentioned such physical processes have a much 
smaller effect on solar neutrino fluxes than the accumu- 



lation of dark matter in the Sun's core. 

In the study of the impact of dark matter in the Sun, 
we will choose to consider an interval of theoretical uncer- 
tainty that takes into account such processes. Therefore, 
considering both the theoretical and experimental uncer- 
tainties, we choose to rule out models that predict a 8 B 
neutrino flux which deviates more than 30% from our 
solar standard model. Similarly, a 7 Be neutrino flux of 
solar models which deviates more than 15% from our so- 
lar standard model can also be ruled out. This is consis- 
tent with the fact that $( 8 B) is two times more sensitive 
to the central temperature than <f>( 7 Be) . The choice of 
this threshold is in agreement with other authors (e.g., 
iTaoso et al.ll2010l) . Furthermore, such intervals of uncer- 
tainties on 8 B and 7 Be also include the solar structure 
variations (e.g. temperature) in the Sun's interior, and 
in particular in the Sun's core that can be attributed to 
some of the physical processes previously mentioned. 

Figure 2 shows the neutrino flux variations in more 
than 30 evolution models of the Sun computed for dif- 
ferent values of rj and (<jv). This result is identi- 
cal for both (s-wave and p-wave) annihilation channels. 
In the r\ vs. (av) plane, the iso-contours correspond- 
ing to $„( 8 £>) and $„( 7 i?e) show neutrino flux varia- 
tions relatively to the solar standard model. We choose 
to present the case of dark matter particles that have 
m x = 10 GeV, spin-dependent scattering cross-section 
a$ d = 10 -40 cm 2 , and spin-independent scattering 
cross-section 175/ = 10~ 36 cm 2 . The neutrino flux varia- 
tion shown follows the decrease of the local halo density 
of particles p x , which follows the relic density fi x specific 
to each set of values (ij, (av)). The observed difference in 
sensitivity between $„( 8 i?) and $> v (' Be) neutrino fluxes, 
is caused by the fact that 8 B neutrinos are produced in 
a more central region than 7 Be neutrinos. 

If we adopt the 30% and 15% fixed thresholds for the 
S B and 7 Be neutrino fluxes, we conclude that we can rule 
out 77— asymmetric dark matter particles with a mass of 
lOGeV, that have an 77 > 10 -12 and an annihilation rate 
(both annihilation channels) larger than 10~ 23 cmV 1 . 
An analysis of other solar neutrino fluxes such as pp neu- 
trinos re-enforces the conclusions reached in this study. 

6. CONCLUSIONS 

In this paper, we proposed a new strategy to study 
the impact of dark matter in the evolution of Sun and 
stars. We started by computing the basic cosmologi- 
cal primitive model that is responsible by the formation 
of the dark matter asymmetric particles, assuming that 
dark matter like baryons has an asymmetry, analogous to 
the baryon asymmetry. As a consequence the population 
of dark matter particles and antiparticles in the current 
Universe is fixed by the measured dark matter density 
fljjMh 2 . The amount of particles and antiparticles de- 
pends on the dark matter asymmetry parameter 77, as 
well as on the dark matter annihilation (p- and s-wave) 
channels. Considering that the Sun is formed in a dark 
matter halo that has the same amount of particles and 
anti-particles, we computed the evolution of the Sun in 
such conditions. 

Using the solar standard model as a reference model 
for the internal structure of the Sun, we have studied the 
impact of asymmetric dark matter on the production of 
solar neutrino fluxes. Following a procedure identical 
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to the analysis proposed in lLopes k, Silkl (|2010al fbT) , we 
describe the impact of asymmetric dark matter on the 
Sun's core. 

Several evolution models of the Sun were computed 
for a range of light dark matter particles. Particles with 
masses smaller than 5 GeV are not considered since evap- 
oration becomes important in this mass range and a large 
number of dark matter particles escape the gravitational 
field of the sta r, significant ly reducing the impact on 
the Sun's core (|Gouldl [1990ft . Similarly, particles with 
a mass above 20 GeV produce a very small dark mat- 
ter core an d their effect in th e Sun's structure is almost 
negligible (|Lopes et al.l 1201 it ). Accordingly, dark mat- 
ter particles with masses between 5 and 20 GeV, pro- 
duce a temperature decrease in the Sun's core that vis- 
ibly affects the neutrino fluxes, although these depend 
on the values of specific parameters, such as the dark 
matter asymmetry and annihilation cross section. The 
impact of more massive dark particles in the evolution of 
a star like the Sun should be significantly smaller than 
in the case discussed in this work. The reason is the 
fact that the concentration of dark matter in the cen- 
tre of the s tar, is usually characterize d by a dark matter 
core radius (|Casanellas fc Lopesj [2009) which is inversely 
proportional to the square root of the mass of the dark 
matter particle. The choice of other parameters of the 
scattering cross section will not change fundamentally 
such results. 

We note incidentally that high annihilation rates for 
light dark matter particles are ruled out if accompa- 
nied by lepton production from observations of Cosmic 
Microwave Background temperature fluctuations and of 
searches for gamma rays from nearby dwarfs. But these 
limits are irrelevant once asymmetric dark matter is im- 
portant, and most importantly for us, assume specific 
annihilation channels. Our solar constraints come from 
accumulation of dark matter particles by the Sun, and 
complement other annihilation limits. The annihilations 
are of course important for fixing the relic abundance. 
However a key point of our paper is that we consider 
asymmetric dark matter for which the dependence on 
specific annihilation channels is very weak. Furthermore, 
the relic abundance is fixed by specifying the annihilation 
rate via freeze-out physics. 

Indeed, freeze-out may be far more complicated than 
given by the simple connection that specifies the "ther- 
mal" freeze-out cross-section in terms of the observed 
relic abundance. Many models violate this condition. 
Our models test independently of freeze-out the very im- 
portant parameter space of annihilation rate today (in 
the low energy universe) whereas relic abundances probe 
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the annihilation rate in a model-dependent way at freeze- 
out (at high energies). Moreover the Sun probes asym- 
metric dark matter at the present solar radius, and com- 
plements the use of dwarfs in the halo and of Cosmic 
Microwave Background in the early universe for testing 
Majorana-type dark matter. 

In this work we have computed the 8 B and 7 Be neu- 
trino flux variations for several 77— parametrised asym- 
metric dark matter particles, corresponding to differ- 
ent values of rj and (crv). The values of 77 and (crv) 
were chosen in such a way that £Idm for each of the 
asymmetric dark matter scenarios considered is consis- 
tent with current observational measurements (cf. Figure 
1). We have presented the results of light dark matter 
particles (typically 10 GeV), with the same character- 
istics as the candidates "observed" by DAM A/LIBRA 
and CoGeNT experiments. We find that such type of 
particles produces large variations on the flux of solar 
neutrinos. The 8 J5 neutrino flux variation takes values 
between 35% and 65% for s-wave and p-wave annihila- 
tion channels, clearly above the uncertainty in the cur- 
rent solar models estimated to be of the order of 15% 
(|Turck-Chieze fc CouvidatJl201ll ). This conclusion holds 
even for our more conservative threshold of 30% on 8 B 
neutrino flux or a threshold of 15% on 7 Be neutrino flux, 
as discussed in the previous section. Likewise, dark mat- 
ter particles with m x = 15 GeV show 8 B neutrino flux 
variations of the order of 16% up to 40%, which can be 
rejected if we consider the 15% threshold of uncertainty 
in the physics of the solar standard model. Therefore, 
it is reasonable to consider that light dark matter parti- 
cles that have a scattering cross section of the order of a 
pico barn, as suggested by several theoretical models to 
explain the DAMA/LIBRA and CoGeNT experiments, 
produce neutrino fluxes that are in disagreement with 
the current neutrino flux observations. On this basis, as- 
suming that stellar evolution is affected by dark matter 
as discussed here, it is reasonable to assume that such 
types of particles cannot exist in the current Universe. 
Forthcoming solar neutrino flux experiments could re- 
strain even further the parameters of the r\— asymmetric 
dark matter candidates proposed by the recent dark mat- 
ter models. 
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